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I Abstract 

The recent LHC indications of a SM-like Higgs boson near 125 GeV are consistent not only 
\ with the Standard Model (SM) but also with Supersymmetry (SUSY). However naturalness 

arguments disfavour the Minimal Supersymmetric Standard Model (MSSM) . We consider the 
Next-to-Minimal Supersymmetric Standard Model (NMSSM) with a SM-like Higgs boson 
near 125 GeV involving relatively light stops and gluinos below 1 TeV in order to satisfy 
naturalness requirements. We are careful to ensure that the chosen values of couplings do 
not become non-perturbative below the grand unification (GUT) scale, although we also 
examine how these limits may be extended by the addition of extra matter to the NMSSM 
at the two-loop level. We then propose four sets of benchmark points corresponding to the 
SM-like Higgs boson being the lightest or the second lightest Higgs state in the NMSSM or 
the NMSSM-with-extra-matter. With the aid of these benchmark points we discuss how the 
NMSSM Higgs boson near 125 GeV may be distinguished from the SM Higgs boson in future 
LHC searches. 
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1 Introduction 



The ATLAS and CMS Collaborations have recently presented the first indication for a Higgs 
boson with a mass in the region ~ 124 — 126 GeV [HE]- An excess of events is observed by the 
ATLAS experiment for a Higgs boson mass hypothesis close to 126 GeV with a maximum local 
statistical significance of 3.6a above the expected SM background and by the CMS experiment 
at 124 GeV with 2.6a maximum local significance. If the ATLAS and CMS signals are combined 
the statistical significance increases, but is still less than the 5a required to claim a discovery. 
Interestingly, the ATLAS signal in the 77 decay channel by itself has a local significance of 
2.8a whereas a SM-like Higgs boson would only have a significance of half this value, leading to 
speculation that the observed Higgs boson is arising from beyond SM physics. In general, these 
results have generated much excitement in the community, and already there are a number of 
papers discussing the implications of such a Higgs boson [3H7]. 

In the Minimal Supersymmetric Standard Model (MSSM) the lightest Higgs boson is lighter 
than about 130-135 GeV, depending on top squark parameters (see e.g. [8] and references 
therein). A 125 GeV SM-like Higgs boson is consistent with the MSSM in the decoupling 
limit. In the limit of decoupling the light Higgs mass is given by 

m 2 h ps Mf cos 2 2/3 + Am 2 h , (1.1) 

where Am 2 is dominated by loops of heavy top quarks and top squarks and tan /3 is the ratio 
of the vacuum expectation values (VEVs) of the two Higgs doublets introduced in the MSSM 
Higgs sector. At large tan /3, we require Am/, ~ 85 GeV which means that a very substantial 
loop contribution, nearly as large as the tree-level mass, is needed to raise the Higgs boson mass 
to 125 GeV. The rather complicated parameter dependence has been studied in [3] where it was 
shown that, with "maximal stop mixing", the lightest stop mass must be mi > 500 GeV (with 
the second stop mass considerably larger) in the MSSM in order to achieve a 125 GeV Higgs 
boson. However one of the motivations for SUSY is to solve the hierarchy or fine-tuning problem 
of the SM [9] . It is well known that such large stop masses typically require a tuning at least of 
order 1% in the MSSM, depending on the parameter choice and the definition of fine-tuning |10j . 

In the light of such fine-tuning considerations, it has been known for some time, even after 
the LEP limit on the Higgs boson mass of 114 GeV, that the fine-tuning of the MSSM could 
be ameliorated in the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [11] . With 
a 125 GeV Higgs boson, this conclusion is greatly strengthened and the NMSSM appears to 
be a much more natural alternative. In the NMSSM, the spectrum of the MSSM is extended 
by one singlet superfield [12H14] (for reviews see [15U16] ). In the NMSSM the supersymmetric 
Higgs mass parameter \x is promoted to a gauge-singlet superfield, S, with a coupling to the 
Higgs doublets, \SH U H^ that is perturbative up to unified scales. In the pure NMSSM values 
of A ~ 0.7 do not spoil the validity of perturbation theory up to the GUT scale only providing 
tan (3 > 4, however the presence of additional extra matter [17] allows smaller values of tan /3 to 
be achieved. The maximum mass of the lightest Higgs boson is 

m 2 h « Aff cos 2 2/3 + AV sin 2 2(3 + Am 2 h (1.2) 

where here we use v = 174 GeV. For \v > Mz, the tree- level contributions to m^ are maximized 
for moderate values of tan f3 rather than by large values of tan /3 as in the MSSM. For example, 
taking A = 0.7 and tan (3 = 2, these tree-level contributions raise the Higgs boson mass to about 
112 GeV, and Am/, > 55 GeV is required. This is to be compared to the MSSM requirement 
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Arrih > 85 GeV. The difference between these two values (numerically about 30 GeV) is signif- 
icant since Ami, depends logarithmically on the stop masses as well as receiving an important 
contribution from stop mixing. This means for example, that, unlike the MSSM, in the case of 
the NMSSM maximal stop mixing is not required to get the Higgs heavy enough. 

The NMSSM has in fact several attractive features as compared to the widely studied MSSM. 
Firstly, the NMSSM naturally solves in an elegant way the so-called fi problem [18 1 of the MSSM, 
namely that the phenomenologically required value for the supersymmetric Higgsino mass [i in 
the vicinity of the electroweak or SUSY breaking scale is not explained. This is automatically 
achieved in the NMSSM, since an effective /i-parameter is dynamically generated when the 
singlet Higgs field acquires a vacuum expectation value of the order of the SUSY breaking scale, 
leading to a fundamental Lagrangian that contains no dimensionful parameters apart from the 
soft SUSY breaking terms. Secondly, as compared to the MSSM, the NMSSM can induce a 
richer phenomenology in the Higgs and neutralino sectors, both in collider and dark matter 
(DM) experiments. For example, heavier Higgs states can decay into lighter ones with sizable 
rates [19H3U] . In addition, a new possibility appears for achieving the correct cosmological relic 
density |31| through the so-called "singlino", i.e. the fifth neutralino of the model, which can 
have weaker-than-usual couplings to SM particles. The NMSSM Higgs boson near 125 GeV may 
also have significantly different properties than the SM Higgs boson 0(6], as we shall discuss 
in some detail here. Thirdly, as already discussed at length above, the NMSSM requires less 
fine-tuning than the MSSM [TT] . 

Since the NMSSM (as in the case for the MSSM) has a rich parameter space, it is convenient 
to consider the standard approach of so-called "benchmark points" in the SUSY parameter space. 
These consist of a few "discrete" parameter configurations of a given SUSY model, which in our 
case are supposed to lead to typical phenomenological features. Using discrete points avoids 
scanning the entire parameter space, focusing instead on representative choices that reflect the 
new interesting features of the model, such as new signals, peculiar mass spectra, and so on. A 
reduced number of points can then be subject to full experimental investigation, without loss of 
substantial theoretical information. While several such benchmark scenarios have been devised 
for the MSSM [32II33] and thoroughly studied in both the collider and the DM contexts, and 
even for the NMSSM |34j . these need to be revised in the light of the 125 GeV Higgs signal. The 
motivation for finding such NMSSM benchmarks is to enable the characteristics of the NMSSM 
Higgs boson near 125 GeV to be identified, so that it may eventually be resolved from that of 
the SM Higgs boson. 

In this paper we present a set of NMSSM benchmark points with a SM-like Higgs boson mass 
near 125 GeV, focussing on the cases where both stop masses are as light as possible in order 
to reduce the fine-tuning, in accordance with the above discussion. The tools to calculate the 
Higgs and SUSY particle spectra in the NMSSM, in particular NMSSMTools, have been available 
for some time [21, 351136] . although these will be supplemented by other codes as discussed later. 
The goal of the paper is to firstly find NMSSM benchmark points which contain a 125 GeV SM- 
like Higgs boson and in addition involve relatively light stops. Since the stops receive radiative 
corrections to their masses from gluinos, we shall also require relatively light gluinos as well 
as having as small an effective /i e g parameter (^ c g = X(S)) as possible [37]. Having relatively 
light stops (and sbottoms), care has to be taken not to be in conflict with direct searches at the 
Tevatron and recent searches at the LHC. In addition, light stops appearing in the triangle loop 
diagrams can significantly affect Higgs production via gluon fusion as well as Higgs decay into 
two photons. Once a parameter set leading to a 125 GeV NMSSM Higgs boson has been found 
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it has to be checked if the production cross-sections and branching ratios are such that they lead 
to a total Higgs production cross-section times branching ratios which is compatible with the 
LHC searches. The obtained branching ratios can be used to ultimately distinguish the NMSSM 
Higgs boson from the SM Higgs state. However, we do not attempt to simulate the LHC search 
strategies, our goal being the much more modest one of providing benchmarks with different 
characteristic types of NMSSM Higgs bosons, which can eventually be studied in more detail. 
The types of benchmark points considered here all involve relatively large values of A > 0.5 
and small values of tan f3 < 3, in order to allow for the least fine-tuning possible involving light 
stops as discussed above. However we are careful to respect the perturbativity bounds on A up 
to the GUT scale, either for the pure NMSSM, or the NMSSM supplemented by three copies 
of extra SU (5) 5 + 5 states, where such bounds are calculated using two-loop renormalisation 
group equations (RGEs). We find that the NMSSM Higgs boson near 125 GeV can come in 
many guises. It may be very SM-like, practically indistinguishable from the SM Higgs boson. 
Or it may have different Higgs production cross-sections and widely varying decay branching 
ratios which enable it to be resolved from the SM Higgs boson. The key distinguishing feature 
of the NMSSM, compared to the MSSM, is the presence of the singlet S which may mix into 
the 125 GeV Higgs mass state to a greater or lesser extent. As the singlet component of the 
125 GeV Higgs boson is increased, either the Hj, or H u components (or both) must be reduced. 
If the Hd component is reduced then this reduces the decay rate into bottom quarks and can 
allow other rare decays like 77 to have larger branching ratios [6]. In addition there is the effect 
of SUSY particle and charged Higgs boson loops in the Higgs coupling to photons which can 
increase or decrease the rate of decays into 77. 

In order to put the present work in context, it is worth to emphasise that this is the first paper 
to propose "natural" (i.e. involving light stops) NMSSM benchmark points with a SM-like Higgs 
boson near 125 GeV. For example, the previous NMSSM benchmark paper [33] was concerned 
with the constrained NMSSM which does not allow light stops consistent with current LHC 
SUSY limits. Also, while this paper was being prepared there have appeared a number of other 
related papers, none of which however are focussed on the task of providing benchmark points. 
For example, the results in [3] were mainly concerned with fine-tuning issues in SUSY models 
with a \SH u Hd coupling for a SM-like Higgs boson near 126 GeV and the actual NMSSM was 
strictly not considered since an explicit fi term was included whereas the trilinear singlet coupling 
was neglected. By contrast another recent paper [5] did consider the actual NMSSM 
although no benchmarks were proposed. As this paper was being finalised, further dedicated 
NMSSM papers with a Higgs boson near 125 GeV have started to appear, in particular [6] 
in which the two photon enhancement was emphasised but without benchmark points, and [7] 
where (versions of) the constrained NMSSM were considered. It should be clear that the present 
paper is both complementary and contemporary to all these papers. 

The layout of the remainder of the paper is as follows. In Section [2] we discuss the MSSM 
and fine-tuning, showing that it leads to constraints on the mass of the heavier stop quark. In 
Section [3] we briefly review the NMSSM, and provide perturbativity limits of the coupling A 
depending on the coupling k as well as tan (3. We show how the limit on A may be increased if 
there is extra matter in the SUSY desert. Section 2] is concerned with constraints from Higgs 
searches, SUSY particle searches and Dark Matter. Section [5] contains the four sets of NMSSM 
Higgs benchmarks near 125 GeV that we are proposing, including the Higgs production cross- 
sections and branching ratios which ultimately will enable it to be distinguished from the SM 
Higgs boson. Section [6] summarises and concludes the paper. Appendix|A] contains the two-loop 



3 



RGEs from which the perturbativity bounds on A were obtained. 



2 The MSSM 

The superpotential of the MSSM is given, in terms of (hatted) superfields, by 
W = fiH u H d + h t Q 3 H u t c R - h b Q 3 H d b c R - KLsHdTR , 



(2.3) 



in which only the third generation fermions have been included (with possible neutrino Yukawa 
couplings have been set to zero), and Q 3 ,L 3 stand for superfields associated with the (t, b) and 
(t,is t ) SU(2) doublets. 

The soft SUSY breaking terms consist of the scalar mass terms for the Higgs and sfermion 
scalar fields which, in terms of the fields corresponding to the complex scalar components of the 
superfields, are given by 



— -Cmass — m H u I H u I + m H d I 

+ ml\Ql\+m}\P R \ + 



- 1 12 I , 2 I t2\ i 2 I ~2 I 

m lR \b R \+m u \L z \+m fR \T R \ , 



and the trilinear interactions between the sfermion and Higgs fields, 

- Aril = BfiH u H d + h t A t Q 3 H u t c R - h b A b Q 3 H d b c R - h T A T L 3 H d f c R + h.c. 
The tree-level MSSM Higgs potential is given by 



(2.4) 



(2.5) 



2 /2 

V = ml\H d \ 2 + ml\H u \ 2 -ml(H d H u + hx.) + ^(H+a a H d + H+a a H u ) 2 + ^-(\H d 



l^l 2 ) 2 
(2.6) 

3/5gi, §2 and g\ are the low energy (GUT normalised) SU(2)w and U(l)y gauge 
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Hu + M and m i 



-Bfi. 



where g' = 

couplings, m\ = rn z Hd + /i 2 , 

In the MSSM, at the 1-loop level, stops contribute to the Higgs boson mass and three 
more parameters become important, the stop soft masses, TOq 3 and mf R , and the stop mixing 
parameter 

X t = A t -ncotp. (2.7) 



The dominant one-loop contribution to the Higgs boson mass depends on the geometric mean 
' J 'Q;"'hr 



of the stop masses, m| = rriQ^m^, and is given by, 
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The Higgs mass is sensitive to the degree of stop mixing through the second term in the brackets, 
and is maximized for \X t \ = X" iax = \f§m~ t , which was referred to as "maximal mixing" above. 

The fine-tuning in the MSSM can be simply understood by examining the leading one-loop 
correction to the Higgs potential, 
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Figure 1: The contribution of one-loop corrections to Eq. H2.11|l as a function of the mass of the heaviest stop m£ 2 
for tan /3 = 10, = 200 GeV, Q = m t = 165 GeV and m h = 300 GeV. Here An = 2-A/M§. Solid, dashed-dotted 
and dashed lines correspond to the mixing angle in the stop sector 9 t = 0,tt/8 and tt/4, respectively. 



where the two stop masses are, 



in 



tl,2 



- ( m% + m~ + 2m 2 t + -M§ cos 2/3 
2 V Q3 tfl 2 



mi - m| + -M^ cos 2/3 - -M§ cos 2/3 ) + 4m 2 X? 



Q3 



(2.10) 



The fine-tuning originates from the fact that AV ^> v 4 , where v = 174 GeV is the combined 
Higgs VEV. By considering the minimization conditions for the Higgs potential, one finds 
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In Eq. (|2.12p Ot is the mixing angle in the stop sector given by 



2m t X t 

sinzffj = — 
(n 

whereas 



sin2fl f = 2 * (2.13) 

v *2 *1' 



Here we set the renormalisation scale Q = mt- From Eq. (|2.11[) one can see that in order to 
avoid tuning, 

A < -M|. (2.14) 

This shows that both stop masses must be light to avoid tuning. For example, defining A// = 
2 • A/M| the absence of any tuning requires A// < 1. This in turn requires the heavier stop 
mass to be below about 500 GeV as illustrated in FigfTJ This constraint on the heavier stop 
mass has not been emphasised in the literature, where often the focus of attention is on the 
lightest stop mass. 

It has been noted that large or maximal stop mixing is associated with large fine-tuning. 
This also follows from FigQ]and Eq. (|2.12|) . Indeed, FigQ] demonstrates that the contribution of 
one-loop corrections to Eq. (|2.1ip increases when the mixing angle in the stop sector becomes 
larger. In fact when B% is close to 7r/4 the last term in Eq. (|2.12|) gives the dominant contribution 
to A enhancing the overall contribution of loop corrections in the minimization condition (|2.11|) 
which determines the mass of the Z-boson. 

Eq. (|2.1ip also indicates that in order to avoid tuning one has to ensure that the parameter 
/i has a reasonably small value. To avoid tuning entirely one should expect [i to be less than 
M-z- However, so small values of the parameter \i are ruled out by chargino searches at LEP. 
Therefore in our analysis we allow the effective fi e Q parameter to be as large as 200 GeV that 
does not result in enormous fine-tuning. 



3 The NMSSM 

In this paper, we only consider the NMSSM with a scale invariant superpotential. Alternative 
models known as the minimal non-minimal super symmetric SM (MNSSM), new minimally- 
extended supersymmetric SM or nearly-minimal supersymmetric SM (nMSSM) or with addi- 
tional U{1)' gauge symmetries have been considered elsewhere [3H], as has the case of explicit 
CP violation [39]. 

The NMSSM superpotential is given, in terms of (hatted) superfields, by 

W = \SH u H d + ^S 3 + h t Q 3 H u t c R - h b Q 3 H d b c R - h T L 3 H d ^ R , (3.15) 

in which only the third generation fermions have been included. The first two terms substitute 
the fiH u H d term in the MSSM superpotential, while the three last terms are the usual gener- 
alization of the Yukawa interactions. The soft SUSY breaking terms consist of the scalar mass 
terms for the Higgs and sfermion scalar fields which, in terms of the fields corresponding to the 
complex scalar components of the superfields, are given by, 

-£ mass = m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 

+ ^IQll+^lill+m^lfell + mljZll + m^lfll . (3.16) 
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The trilinear interactions between the sfermion and Higgs fields are, 



- £ tri i = XA x H u H d S + ^kA k S 3 + h t A t Q 3 Hj c R - h b A b Q 3 H d b c R - h T A T L 3 H d f R + h.c. .(3.17) 

In the unconstrained NMSSM considered here, with non-universal soft terms at the GUT scale, 
the three SUSY breaking masses squared for H u , H d and S appearing in -Anass can be expressed in 
terms of their VEVs through the three minimization conditions of the scalar potential. Thus, in 
contrast to the MSSM (where one has only two free parameters at the tree level, generally chosen 
to be the ratio of Higgs vacuum expectation values, tan /3, and the mass of the pseudoscalar Higgs 
boson), the Higgs sector of the NMSSM is described by the six parameters 

A , k , A x , A K , tan/3 = (H u )/(H d ) and /x cfr = \{S) . (3.18) 

We follow the sign conventions such that the parameters A and tan (3 are positive, while the 
parameters k, A\, A k and /i e ff can have both signs. 

In addition to these six parameters of the Higgs sector, one needs to specify the soft SUSY 
breaking mass terms in Eq. (|3.16p for the scalars, the trilinear couplings in Eq. (|3.17p as well as 
the gaugino soft SUSY breaking mass parameters to describe the model completely, 



C - - 

'—gauginos — 2 



MxBB + M 2 ^W a W a + M 3 J2G a G a + h.c. 



a=l a=l 



(3.19) 



Clearly, in the limit A — > with finite (i e s, the NMSSM turns into the MSSM with a decoupled 
singlet sector. Whereas the phenomenology of the NMSSM for A — > could still differ somewhat 
from the MSSM in the case where the lightest SUSY particle is the singlino (and hence with the 
possibility of a long lived next-to- lightest SUSY particle [ID] ) , we will not consider this situation 
here. In fact we shall be interested exclusively in large values of A (i.e. A > 0.5) in order to 
increase the tree- level Higgs mass as in Eq. (|1.2p . For the same reason we shall also focus on 
moderate values of tan /3 (tan ft = 2, 3) that result in the relatively large values of the top quark 
Yukawa coupling ht at low energies. 








0.1 


0.2 


0.3 


0.4 


0.5 


tan = 2 


0.62 


0.61 


0.60 


0.58 


0.53 


0.42 


tan /3 = 3 


0.68 


0.68 


0.66 


0.63 


0.56 


0.45 



Table 1: Two-loop upper bounds on A(Afz) for different values of n(Mz) and tan/3 in the NMSSM. 



The growth of the Yukawa couplings ht, A and k at the electroweak (EW) scale entails the 
increase of their values at the Grand Unification scale Mj resulting in the appearance of the 
Landau pole. Large values of h t (Mx), X(Mx) and k(Mx) spoil the applicability of perturbation 
theory at high energies so that the RGEs cannot be used for an adequate description of the 
evolution of gauge and Yukawa couplings at high scales Q ~ Mx- The requirement of validity 
of perturbation theory up to the Grand Unification scale restricts the interval of variations 
of Yukawa couplings at the EW scale. In particular, the assumption that perturbative physics 
continues up to the scale Mx sets an upper limit on the low energy value of \{Mz) for each fixed 
set of k(Mz) and ht{M t ) (or tan/3). With decreasing (increasing) k(Mz) the maximal possible 
value of X(Mz), which is consistent with perturbative gauge coupling unification, increases 
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(decreases) for each particular value of tan j3. In Table Q] we display two-loop upper bounds 
on X(Mz) for different values of k(Mz) and tan/3 in the NMSSM. As one can see the allowed 
range for the Yukawa couplings varies when tan/3 changes. Indeed, for tan/3 = 2 the value of 
X(Mz) should be smaller than 0.62 to ensure the validity of perturbation theory up to the scale 
Mx- At large tan/3 the allowed range for the Yukawa couplings enlarges. The upper bound on 
\(Mz) grows with increasing tan (3 because the top-quark Yukawa coupling decreases. At large 
tan (3 (i.e. tan j3 > 4 — 5) the upper bound on X(Mz) approaches the saturation limit where 
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Table 2: Two-loop upper bounds on \{Mz) for different values of k(A1z) and tan/3 within the NMSSM supple- 
mented by three SU (5) (5 + 5)-plets. The masses of all extra exotic states are set to be equal to 300 GeV. 
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Table 3: Two-loop upper bounds on X(Mz) for different values of k(Mz) and tan/3 within the NMSSM supple- 
mented by three SU (5) (5 + 5)-plets. The masses of all extra exotic states are set to be equal to 1 TeV. 



The renormalisation group (RG) flow of the Yukawa couplings depends rather strongly on 
the values of the gauge couplings at the intermediate scales. To demonstrate this we examine 
the RG flow of gauge and Yukawa couplings within the SUSY model that contains three extra 
577(5) (5 + 5)-plets that survive to low energies and can form three 10-plets in the SUSY-GUT 
model based on the 50(10) gauge group. In this SUSY model the strong gauge coupling has 
a zero one-loop beta function whereas at two-loop level the coupling has a mild growth as the 
renormalisation scale increases. Since extra states form complete SU(5) multiplets the high- 
energy scale where the unification of the gauge couplings takes place remains almost the same 
as in the MSSM. At the same time extra multiplets of matter change the running of the gauge 
couplings so that their values at the intermediate scale rise substantially. In fact the two-loop 
beta functions of the SM gauge couplings are quite close to their saturation limits when these 
couplings blow up at the GUT scale. Further enlargement of the particle content can lead to 
the appearance of the Landau pole during the evolution of the gauge couplings from Mz to 
Mx- Because the beta functions are so close to the saturation limits the RG flow of the gauge 
couplings (i.e. their values at the intermediate scale) also depends on the masses of the extra 
exotic states. Since gi{Q) occurs in the right-hand side of the RGEs for the Yukawa couplings 
with negative sign the growth of the gauge couplings prevents the appearance of the Landau 
pole in the evolution of these couplings. It means that for each value of ht(Mt) (or tan/3) and 
k(Mz) the upper limit on X(Mz) increases as compared with the NMSSM. The two-loop upper 
bounds on X{Mz) for different values of k(Mz) and tan /3 in the NMSSM supplemented by three 
5O(10) 10-plets, or equivalently three SU(5) (5 + 5)-plets are displayed in Table [2] for the case 
that the masses of all extra exotic states are set to be equal to 300 GeV and in Table [3] in case 
they are set to be equal to 1 TeV. The two- loop RGEs used to obtain these results are given 
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in Appendix E Because the RG flow of the gauge couplings depends on the masses of extra 
exotic states the upper bounds on \(Mz) presented in Tables [2] and are slightly different. 
The restrictions on \{Mz) obtained in this section are useful for the phenomenological analysis 
which we are going to consider next. 

4 Constraints from Higgs Boson Searches, SUSY Particle Sear- 
ches and Dark Matter 

Our scenarios are subject to constraints on the Higgs boson masses from the direct searches at 
LEP, Tevatron and the LHC. Also the SUSY particle masses have to be compatible with the 
limits given by the experiments. Finally, the currently measured value of the relic density shall 
be reproduced. Further constraints arise from the low-energy observables. 



4.1 Higgs boson searches 

We start by discussing the constraints which arise from the LHC search for the Higgs boson. 
At the LHC, the most relevant Higgs boson production channels for neutral (N)MSSM Higgs 
bosons are given by gluon fusion, gauge boson fusion, Higgs-strahlung and associated production 
with a heavy quark pair. The two main mechanisms for charged Higgs boson production are top 
quark decay and associated production with a heavy quark pair. For reviews, see [81. 1411. 1^2] . As 
in our scenarios the charged Higgs boson mass is larger than 450 GeV and hence well beyond the 
sensitivity of Tevatron and current LHC searches, we will discuss in the following only neutral 
Higgs boson production. 

Gluon fusion In the SM and in SUSY extensions, such as the (N)MSSM, for low values of 
tan/3, the most important production channel is given by gluon fusion [43] . In the NMSSM we 
have 



gg -> H % and gg -»■ A s , i = 1, 2, 3, j = 1, 2 . (4.20) 

Since this is the dominant Higgs production mechanism for a 125 GeV Higgs boson at the LHC, 
we find it convenient to define for later use the ratio of the gluon fusion production cross-section 
for the Higgs boson Hi in the NMSSM to the gluon fusion production cross-section for a SM 
Higgs boson H SkI of same mass as Hi, 

Gluon fusion is mediated by heavy quark loops in the SM and additionally by heavy squark loops 
in the (N)MSSM. It is subject to important higher-order QCD corrections. For the SM, they 
have been calculated at next-to-leading order (NLO) [44] including the full mass dependence 
of the loop particles and in the heavy top quark limit, and up to next-to-next-to-leading order 
(NNLO) in the heavy top quark limit [15]. The cross-section has been improved by soft-gluon 
resummation at next-to-next-to- leading logarithmic (NNLL) accuracy [16]. Top quark mass 
effects on the NNLO loop corrections have been studied in [47] . and the EW corrections have 
been provided in [48J. In the MSSM, the QCD corrections have been calculated up to NLO [44] . 
The QCD corrections to squark loops have been first considered in [19] and at full NLO SUSY- 
QCD in the heavy mass limit in [SD]. The (s)bottom quark contributions at NLO SUSY-QCD 
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have been taken into account through an asymptotic expansion in the SUSY particle masses [51] . 
For squark masses below ~ 400 GeV, mass effects play a role and can alter the cross-section by 
up to 15% compared to the heavy mass limit as has been shown for the QCD corrections to the 
squark loops in [521153] - The SUSY QCD corrections including the full mass dependence of all 
loop particles have been provided by [S3]- The mass effects turn out to be sizeable. The NNLO 
SUSY-QCD corrections from the (s)top quark sector to the matching coefficient determining the 
effective Higgs gluon vertex have been calculated in [55J. 

The gluon fusion cross-section has been implemented in the Fortran code HIGLU |56| up 
to NNLO QCD. While at NLO the full mass dependence of the loop particles is taken into 
account the NNLO corrections are obtained in an effective theory approach. In the MSSM the 
full squark mass dependence in the NLO QCD corrections to the squark loops is included [53] . 
Note, however, that in the MSSM at NNLO the mismatch in the QCD corrections to the effective 
vertex is not taken into account, neither the SUSY QCD corrections to the effective vertex. The 
former should be only a minor effect, though, as the dominant effect of the QCD higher-order 
corrections stems from the gluon radiation. Furthermore, the EW corrections to the SM can be 
obtained with HIGLU. In order to check if our scenario is compatible with the recent LHC results, 
we need the cross-section of a SM-like Higgs boson of 124 to 126 GeV. The experiments include 
in their analyses the NNLO QCD (CMS also the NNLL QCD and NLO EW) corrections [HE] to 
the gluon fusion cross-section as provided by the Higgs Cross-Section Working Group [42] . For 
SUSY, however, the EW corrections are not available. In order to be consistent, we therefore 
compare in the following the NMSSM cross-section to the SM cross-section at NNLO QCD. As 
the QCD corrections are not affected by modifications of the Higgs couplings to the (s)quarks, 
the NMSSM cross-section can be obtained with the program HIGLU by multiplying the MSSM 
Higgs couplings with the corresponding modification factor of the NMSSM Higgs couplings to 
the (s)quarks with respect to the MSSM case. We have implemented these coupling modification 
factors in the most recent HIGLU version 3.11. 

W/Z-hoson fusion Gauge boson fusion [57J plays an important role for light CP-even Higgs 
boson production in "the SM limit Q 

qq->qq+W*W*/Z*Z*->qqHi, i = 1,2,3. (4.22) 

Otherwise the (N)MSSM cross-section is suppressed with respect to the SM case by mixing 
angles entering the Higgs couplings to the gauge bosons. The NLO QCD corrections are of 
0(10%) of the total cross-section [HH1I53] . The full EW and QCD corrections to the SM are 
0(5%) [60]. The NNLO QCD effects on the cross-section amount to ~ 2% [51]. The SUSY 
QCD and SUSY EW corrections are small [62,63]. Once again, as QCD corrections are not 
affected by the Higgs couplings to the gauge bosons, the QCD corrected NMSSM gauge boson 
fusion production cross-sections can be derived from the QCD corrected SM cross-section by 
simply applying the modification factor of the respective NMSSM Higgs coupling to the gauge 
bosons with respect to the SM coupling, 

^ D SSM (qqH t ) = (j™*.) 2 a™ D (qqH™) , V = W, Z , (4.23) 
\9VVHSM J 

where g denotes the coupling. The EW corrections, however, cannot be taken over. We have 
obtained the SM production cross-section at NLO QCD from the program VV2H [64]. While the 



lr The quark q stands for a generic quark flavour, which is different for the two quarks in case of W^-boson 
fusion. The same notation is applied below in Higgs-strahlung. 
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experiments use the SM cross-section at NNLO QCD (CMS also at NLO EW), the effects of 
these additional corrections in the SM limit, where we compare our NMSSM Higgs cross-section 
to the SM case, are small. 

Higgs-strahlung The CP-even Higgs bosons can also be produced in Higgs-strahlung [65:, 

qq^VHi, V = W,Z, 1 = 1,2,3, (4.24) 

with the NMSSM cross-section always being suppressed by mixing angles compared to the SM 
cross-section. The QCD corrections apply both to the SM and (N)MSSM case. While the NLO 
QCD corrections increase the cross-section by 0(30%) [59 , 66J the NNLO QCD corrections are 
small [67j- The full EW corrections are only known for the SM and decrease the cross-section by 
0(5 - 10%) [68]. The SUSY-QCD corrections amount to less than a few percent [62]. The NLO 
QCD SM Higgs-strahlung cross-section has been obtained with the program V2HV [M]. The 
NMSSM Higgs production cross-sections can be derived from it by applying the Higgs coupling 
modification factors, 

a^^{V Hi ) = (9vmL\ 2 a s Q UVH SM ) , V = W,Z. (4.25) 
\9VVH SM J 

The experiments use the QCD corrected cross-section up to NNLO (CMS also including the 
NLO EW corrections). While we neglect the NNLO and EW corrections, we do not expect this 
to influence significantly the total cross-section composed of all production channels, in view of 
the small size of the total Higgs-strahlung cross-section itself. 

Associated production with heavy quarks Associated production of (N)MSSM Higgs bo- 
sons with top quarks [69j only plays a role for the light scalar Higgs particle and small values of 
tan/3 due to the suppression of the Higgs couplings to top quarks ~ 1/ tan/3. While associated 
production with bottom quarks [69 ,70 does not play a role in the SM, in the (N)MSSM this cross- 
section becomes important for large values of tan /3 and can exceed the gluon fusion cross-section. 
As our scenarios include small values of tan (3 we will not further discuss this cross-section here. 
The values for the SM ttH cross section including NLO QCD corrections [71], which are of 
moderate size, can be obtained from the Higgs Cross Section Working Group homepage [72] . 
From these we derived the NMSSM cross-section values by replacing the SM Yukawa couplings 
with the NMSSM Yukawa couplings, 

a^ D SSM (tt Hl ) = (J**.) 2 *$ D (tiH SM ) ■ (4.26) 

\9ttHSM J 

The NLO SUSY QCD corrections, which have not been taken into account by the experiments, 
are of moderate size [73] . 



4.2 Constraints from the LHC Searches 

Recent results presented by the ATLAS [1] and the CMS [2] Collaborations seem to indicate 
a Higgs boson of mass of 126 and 124 GeV, respectively. Based on the dataset corresponding 
to an integrated luminosity of up to 4.9 fb _1 collected at yfs = 7 TeV, an excess of events 
is observed by the ATLAS experiment for a Higgs boson mass hypothesis close to Mu = 126 
GeV with a maximum local significance of 3.6a above the expected SM background. The three 
most sensitive channels in this mass range are given by H — > 77, H — > ZZ^ — > l + l~l + l~ and 
H —7- WW^*) —> l + vl~v. The CMS Collaboration presented results of SM Higgs boson searches 
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in the mass range 100-600 GeV in 5 decay modes, H — >■ 77, bb, tt, WW and ZZ. The data 
correspond to an integrated luminosity of up to 4.7 fb _1 at yfs = 7 TeV. A modest excess of 
events is observed for Higgs boson mass hypotheses towards the low end of the investigated Higgs 
mass range. The maximum local significance amounts to 2.6a for a Higgs boson mass hypothesis 
of Mu = 124 GeV. For our NMSSM benchmark scenarios presented below to be consistent with 
these LHC results we demand the production cross-section of the SM-like NMSSM Higgs boson 
with mass 124 GeV to 126 GeV (depending on the scenario) to be compatible within 20% 
with the production cross-section of a SM Higgs boson of same mass. The 20% are driven by 
the theoretical uncertainty on the inclusive Higgs production cross-section given by the sum 
of the most relevant production channels at low values of tan /3, i.e. gluon fusion, weak boson 
fusion, Higgs-strahlung and tt Higgs production. The theoretical error is largest for the gluon 
fusion cross-section with 10-15% at these Higgs mass values and yfs = 7 TeV [32], and which 
contributes dominantly to the inclusive production. We do not consider any experimental error 
since this is beyond our scope. 

For simplicity, and since these search channels are common to both experimental analyses, 
we consider the Higgs decays into 77, ZZ — > Al and WW — > 2l2u. In order to get an estimate 
of how closely the NMSSM Higgs resembles the SM Higgs in LHC searches in these channels we 
define the ratios of branching ratios into massive gauge boson final states VV, where V = W, Z$ 
and into 77, respectively, for an NMSSM Higgs boson Hi and the SM Higgs boson H SM of same 
mass, 

BRjH^VV) BR{Hj — > 77) 

Rwm = BR{HSM - yv) and = br{hSM ^ ^ . (4.27) 

We also define analogously Rr tot for the total widths, 

^^"fSP^y' (4 - 28) 

and R b ^(Hi) for the decay into bb. Although this final state is not useful for LHC searches, 
it is interesting to show as in this mass range and for small values of tan f3 the decay into bb 
contributes dominantly to the total width. Depending on of how much Hd component is in the 
mass eigenstate Hi it is enhanced or suppressed compared to the SM. This directly influences 
the total width and hence the branching ratios of the other final states. 

For a crude estimate of the total Higgs cross-section at the LHC, we can combine these 
channels in quadrature in contrast to the experiments which do a sophisticated statistical com- 
bination of the various search channels, which is, however, beyond the scope of our theoretical 
analysis. Our results should therefore only be regarded as a rough estimate which is indicative 
enough, however, at the present status of the experimental research, to exclude or not exclude a 
benchmark scenario. We hence demand for a scenario to be valid that one of the NMSSM Higgs 
bosons Hi satisfies 

0.8a tot (H SM ) < a tot (Hi) < 1.2a tot (H SM ) , (4.29) 

where 

a tot {H) = a ind (H)^BR 2 (H ^ 77) + 16 5i? 2 (F ZZ)BR 2 (Z -> ll)BR 2 (Z -> 11) 

+ WBR 2 (H ^WW)BR 2 (W ^lv)BR 2 (W ^lv)} 1/2 . (4.30) 



2 The ratio is the same for WW and ZZ final states, respectively, as the NMSSM coupling to WW and ZZ is 
suppressed by the same factor compared to the SM. 
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The inclusive cross-section <Ji nc i is composed of gluon fusion, vector boson fusion, Higgs-strahlung 
and associated production with tt, 

a incl {H) = a{gg -> H) + a(Hqq) + a(WH) + a{ZH) + a(tiH) , (4.31) 

with if = Hi,H SM , respectively, subject to the constraint M h sm = Mn t = 124-126 GeV, 
depending on the scenario under consideration. It is dominated by the gluon fusion cross- 
section. The factors 16 in Eq. (|4,30p arise from the sum of the four possible lepton final states 
in the decays of the Z- and FU-boson pairs, respectively. (We neglect interference effects.) For 
the gauge boson branching ratios we chose the values given by the Particle Data Group |74j . 

BR(Z 11) = 0.0366 , BR{W Iv) = 0.1080 . (4.32) 

It is useful to define 

R<7 tot (Hi) = at *j H s M y ( 4 - 33 ) 

in order to provide a measure of how closely the NMSSM Higgs resembles the SM Higgs in 
the most important current LHC search channels. The total cross-section is dominated by the 
VUVF-boson final state due to the large branching ratio. As we will see below, in the NMSSM 
the branching ratio into 77 can be enhanced for certain parameter configurations compared to 
the SM. To illustrate this effect, we therefore also calculate separately the ratios of the expected 
number of events in the NMSSM compared to the SM for the 77 final state and for the VV 
(V = W, Z) final state, which is the same for V = W or Z . They are given by 

R amc[ (H t )R^(H t ) and R aincl (Hi)Rvv{Hi) , (4.34) 

where R aincl {H t ) = <j inc i{Hi)/<j ind {H SM ). 



4.3 NMSSM spectrum and NMSSM Higgs boson branching ratios 

The SUSY particle and NMSSM Higgs boson masses and branching ratios are calculated with 
the program package NMSSMTools [211 [35l [36] . As for the NMSSM Higgs boson masses, the 
leading one- loop contributions due to heavy (s) quark loops calculated in the effective potential 
approach [75], the one-loop contributions due to chargino, neutralino and scalar loops in leading 
logarithmic order in Ref. [76j and the leading logarithmic two-loop terms of 0(ata s ) and O(af), 
taken over from the MSSM results, have been implemented in NMSSMTools. The full one- 
loop contributions have been computed in the DR renormalisation scheme |77^ 78j and also in 
a mixed on-shell (OS) and DR scheme as well as in a pure OS scheme [79] . Furthermore, the 
0(ata s + a.ba s ) corrections have been provided in the approximation of zero external momentum 
[77] . The corrections provided by Ref. [77] have been implemented in NMSSMTools as well. 

The calculation of the NMSSM Higgs boson decay widths and branching ratios within 
NMSSMTools is performed by the Fortran code NMHDECAY [2HES] which uses to some extent 
parts of the Fortran code HDECAY [80l[81] that calculates SM and MSSM Higgs boson partial 
widths and branching ratios. The calculation of the SUSY particle branching ratios on the 
other hand with the Fortran code NMSDECAY [82] is based on a generalisation of the Fortran code 
SDECAY [8TII83] to the NMSSM case. NMSSMTools provides the output for the complete NMSSM 
particle spectrum and mixing angles and for the decays in the SUSY Les Houches format |84j . 
The latter can be read in by our own Fortran version for NMSSM Higgs boson decays based on 
an extension of the latest HDECAY version. It reads in the particle spectrum and mixing angles 
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created with NMSSMTools, calculates internally the NMSSM Higgs boson couplings and uses 
them to calculate the Higgs decay widths and branching ratios. The results for the branching 
ratios from NMSSMTools and from our own program agree reasonably well and the differences 
in the total cross-section Eq. (|4.30p . obtained with the results from the two programs, due to 
deviations in the branching ratios are in the percent range. 

4.4 Constraints from Dark Matter, Low Energy Observables, LEP and Teva- 



Based on an interface between NMHDECAY and micrOMEGAs [8SJ the relic abundance of the NMSSM 
dark matter candidate Xi can be evaluated using NMSSMTools. As an independent check, we 
also used the stand alone code micrOMEGAs to calculate the relic density. All the relevant 
cross-sections for the lightest neutralino annihilation and co-annihilation are computed. The 
density evolution equation is numerically solved and the relic density of Xi 1S calculated. The 
differences in the result for the relic density calculated with both tools are negligible. The results 
are compared with the "WMAP" constraint 0.094 < Qcdm/i 2 < 0.136 at the 2a level [86] . 




When the spectrum and the couplings of the Higgs and SUSY particles are computed with 
NMSSMTools, constraints from low-energy observables as well as available Tevatron and LEP 
constraints are checked. The results of the four LEP collaborations, combined by the LEP Higgs 
Working Group, are included [87] • More specifically, the following experimental constraints are 
taken into account: 

(i) The masses of the neutralinos as well as their couplings to the Z° boson are compared 
with the LEP constraints from direct searches and from the invisible Z° boson width. 

(ii) Direct bounds from LEP and Tevatron on the masses of the charged particles (h^, x^i 
q, I) and the gluino are taken into account. 

(iii) Constraints on the Higgs production rates from all channels studied at LEP. These include 
in particular ZHi production, Hi being any of the CP-even Higgs particles, with all 
possible two body decay modes of Hi (into b quarks, r leptons, jets, photons or invisible), 
and all possible decay modes of Hi of the form Hi —> AjAj, Aj being any of the CP-odd 
Higgs particles, with all possible combinations of Aj decays into b quarks, c quarks, r 
leptons and jets. Also considered is the associated production mode e + e~ — > HiAj with, 
possibly, Hi — > AjAj. (In practice, for our purposes, only combinations of i = 1,2 with 
j = 1 are phenomenologically relevant.) 

(iv) Experimental constraints from B physics [88] such as the branching ratios of the rare 
decays BR(i3 — > Xtf), BR(B S — > /^ + /U~) and BR(£> + — > t + u t ) and the mass differences 
AM S and AM^, are also implemented; compatibility of each point in parameter space with 
the current experimental bounds is required at the two sigma level. 

4.5 Constraints on SUSY particle masses from the LHC 

The ATLAS [89] and CMS [90] searches in final states with jets and missing transverse energy 
EJp' lss , with large jet multiplicities and E™ tss , with heavy flavour jets and EJp iss within simplified 
models and mSUGRA/constrained MSSM (CMSSM) models set limits on the masses of gluino 
and squark masses. Further constraints are obtained from searches in final states with leptons 
and taus [91J. The precise exclusion limits depend on the investigated final state, the value of 




the neutralino and/or chargino masses and the considered model. Light gluino (below about 
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600 GeV) and squark masses (below about 700 GeV) are excluded. The limits cannot be 
applied, however, to the third generation squarks in a model-independent way. Recent analyses 
scanning over the physical stop and sbottom masses and translating the limits to the third 
generation squark sector have shown that the sbottom and stop masses can still be as light as 
~ 200 — 300 GeV depending on the details of the spectrum [92J. Especially scenarios, where 
the lightest stop is the next-to-lightest SUSY particle (NLSP) and nearly degenerate with the 
lightest neutralino assumed to be the lightest SUSY particle (LSP), are challenging for the 
experiments. Such scenarios can be consistent with Dark Matter constraints due to possible 
co-annihilation [93]. If the t\ — X\ mass difference is small enough, the flavour-changing neutral 
current decay t\ — > ex® [93] is dominating and can compete with the four-body decay into the 
LSP, a b quark and a fermion pair |95j . Limits have been placed by the Tevatron searches [SB] 
and depending on the neutralino mass still allow for very light stops down to 100 GeV. The 
authors of Ref. [97| found that translating the LHC limits to stop searches in the co-annihilation 
scenario |98_H allows for stops lighter than ~ 400 GeV down to 160 GeV. We are not aware, 
however, of any dedicated LHC analysis which excludes very light stops. 

5 The benchmark points 

In this section we present benchmark points for the NMSSM with a SM-like Higgs boson near 125 
GeV. The Higgs sector of the NMSSM has a rich parameter space including A, k, A\, A k , tan (3 
and the effective [i parameter. According to the SLHA format [84] . these parameters are un- 
derstood as running DR parameters taken at the SUSY scale M = 1 TeV while tan /3 is taken 
at the scale of the Z boson mass, Mz- In order not to violate tree-level naturalness we set 
Me// — 200 GeV for all the considered points. The Higgs sector is strongly influenced by the 
stop sector via radiative corrections where we further need to specify the soft SUSY breaking 
masses Mq^ , M^ R and the mixing parameter Xt defined in Eq. (|2.7p . The main advantage of 
the NMSSM over the MSSM, regarding a SM-like Higgs boson near 125 GeV, is that the stop 
masses are allowed to be much lighter, making the NMSSM much more technically natural than 
the MSSM. Thus all of the benchmarks discussed here will involve relatively light stops, with 
masses well below 1 TeV. We choose low values of tan f3 = 2, 3 in order to maximise the tree-level 
contribution to the Higgs boson mass, allowing the stops to be lighter. For very light stops, in 
order not to be in conflict with the present exclusion limits, the difference between the t\ and 
Xi masses should be less than ~ 20 — 30 GeV to be in the co-annihilation region. By choosing 
the right-handed stop to be the lightest top squark the sbottoms are still heavy enough to fulfill 
the LHC limits of about 300 GeV |100] also in this case. On the other hand, the remaining 
squarks and sleptons may be heavier without affecting fine-tuning. In order to satisfy in partic- 
ular the LHC search limits for the squarks of the first two generations, we set all their masses 
to be close to 1 TeV and, for simplicity, also those of the sleptons. To be precise, for the first 
and second squark and slepton families and for the stau sector we always take the soft SUSY 
breaking masses and trilinear couplings to be 1 TeV, and furthermore, the right-handed soft 
SUSY breaking mass and trilinear coupling of the sbottom sector is set to 1 TeV, i.e. (U = u, c, 
D = d,s, E = e, fj,, t) 



A v = A D = A b = A E = lTeV. 




= M& =M f 



= 1 TeV, 



(5.35) 



3 For stop searches in scenarios with light gravitinos see |99| . 
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This results in physical masses of ~ 1 TeV for the first and second family squarks and sleptons 
as well as the heavier sbottom. These masses can readily be increased without affecting the 
properties of the quoted benchmark points appreciably. The gaugino mass parameters have 
been set such that they fulfill roughly GUT relations. Special attention has been paid, however, 
not to choose the gluino mass too heavy in order to avoid fine-tuning. Before discussing the 
benchmark points, a few technical preliminaries are in order. The masses for the Higgs bosons 
and SUSY particles have been obtained using NMSSMTools. In the presentation of our benchmark 
scenarios, for the SM-like Higgs boson we furthermore include the ratios of the branching ratios, 
c/. Eq. (023), into 77, bb and VV (V = Z, W) as well as the ratio of the total widths Eq. flOED 
for the NMSSM and SM Higgs boson having the same mass. The NMSSM branching ratios 
and total width have been obtained with NMSSMTools and cross-checked against a private code 
whereas the SM values have been calculated with HDECAY. The latter are shown separately in 
Table 0] for a SM Higgs boson with the mass values corresponding to the various masses of the 
SM-like NMSSM Higgs boson in our benchmark scenarios. The NMSSM values can be obtained 
by multiplication with the corresponding ratios presented in the benchmark tables, although it 
is mainly their relative values, compared to the SM, that concern us here. Note that HDECAY 
includes the double off-shell decays into massive vector bosons, whereas NMSSMTools does not. 
We therefore turned off the double off-shell decays also in HDECAY. This explains why the SM 
values given in Table 2] differ from the values given on the website of the Higgs Cross Section 
Working Group |72) . There are further differences between the two programs in the calculation 
of the various partial widths. Thus HDECAY includes the full NNNLO corrections to the top loops 
in the decay into gluons. Also it uses slightly different running bottom and charm quark masses. 
Taking all these effects into account we estimate the theoretical error on the ratios of branching 
ratios, which are calculated with these two different programs, to be of the order of 5%. This 
should be kept in mind when discussing the benchmark scenarios. 

Table 4: Branching ratios into 77, ZZ, WW , bb and total width of a SM Higgs boson of mass between 123.5 and 
126.5 GeV. 



M„ [GeV] 


BR(H -> 77) 


BR(H ->■ ZZ) 


BR(# -> WW) 


BR(# -)■ bb) 


Ttot [GeV] 


123.5 


2.334 • 1(T 3 


0.018 


0.174 


0.616 


3.773 • 10 -3 


123.6 


2.335 • 10" 3 


0.018 


0.175 


0.615 


3.784 • 10 -3 


123.7 


2.336 • 1(T 3 


0.019 


0.177 


0.613 


3.795 • 10" 3 


123.8 


2.337- 1(T 3 


0.019 


0.179 


0.612 


3.807 • 10" 3 


124 


2.338 ■ 1(T 3 


0.019 


0.182 


0.609 


3.830 • 10" 3 


124.5 


2.342 ■ 10" 3 


0.020 


0.189 


0.601 


3.890 • 10" 3 


124.6 


2.343 • 10" 3 


0.021 


0.191 


0.600 


3.903 • 10" 3 


125 


2.345 • 1(T 3 


0.021 


0.197 


0.594 


3.953 • 10~ 3 


125.8 


2.347 • 1(T 3 


0.023 


0.210 


0.581 


4.058 • 10~ 3 


126 


2.348 • 1(T 3 


0.024 


0.213 


0.578 


4.085 • 10~ 3 


126.2 


2.348 • 1(T 3 


0.024 


0.217 


0.575 


4.113 • 10~ 3 


126.5 


2.348 • 1(T 3 


0.025 


0.222 


0.570 


4.155 • 10~ 3 



As mentioned above, in order to be compatible with the recent LHC results for the Higgs 
boson search we demand the total cross-section as defined in Eq. (|4.30p to be within 20 % equal 
to the corresponding SM cross-section. We therefore include in our tables for the benchmark 
points the ratio R atot , Eq. (|4.33[) . of the NMSSM and SM total cross-section for = 7 TeV 
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and for completeness also the ratio -R CTfl9 , Eq. ()4.2ip . of the NMSSM and SM gluon fusion cross- 
sections since gluon fusion is the dominant contribution to inclusive Higgs production for low 
values of tan j3 at the LHC. The latter has been obtained with the Fortran code HIGLU at NNLO 
QCD and includes the squark loop contributions. Note, that we did not include electroweak 
corrections. Furthermore, we explicitly verified that R atot is practically the same using NLO 
or NNLO QCD gluon fusion cross-sections. The values of the gluon fusion cross-section and of 
the total cross-section at NLO and NNLO QCD are shown separately for the SM in Table [5) 
However, the ratios for the cross-sections which we present in the tables, Ra gg , R<j to ti are f° r the 
case of gluon fusion production at NNLO QCD. With the total cross-section being dominated 
by the VFVF-boson final state, in order to illustrate interesting effects in the NMSSM branching 
ratios compared to the SM ones, for the SM-like Higgs boson we also give separately the ratios of 
the expected number of events in the 77 and VV (V = W, Z) final states as given by Eq. (|4.34p . 

All our presented scenarios fulfill the cross-section constraint Eq. (|4.29p and are hence com- 
patible with the LHC searches, keeping in mind though that as theorists we can do only a rough 
estimate here. Furthermore, they fulfill the constraints from low-energy parameters as specified 
in section 0] and are compatible with the measurement of the relic density. In principle they 
could also account for the 3a deviation of the muon anomalous magnetic momentum from the 
SM if we were to lower the smuon mass, but for clarity we have taken all first and second family 
squark and slepton masses to be close to 1 TeV, as discussed above. 

We consider four different sets of NMSSM benchmark points as follows: 

• NMSSM with Lightest Higgs being SM-like near 125 GeV 

This is achieved with A = 0.57 — 0.64 and k = 0.18 — 0.25 such that A does not blow up 
below the GUT scale in the usual NMSSM with no extra matter as in Table [TJ This set 
of benchmarks is displayed in Table [H 

• NMSSM with Second Lightest Higgs being SM-like near 125 GeV 

This is achieved with A = 0.55 — 0.67 and k = 0.10 — 0.31 such that A does not blow up 
below the GUT scale in the usual NMSSM with no extra matter as in Table [TJ This set 
of benchmarks is displayed in Table \7\ 

• NMSSM-with-extra-matter and Second Higgs being SM-like near 125 GeV 

This is achieved with A = 0.68 — 0.69 and k = 0.06 — 0.20 such that A does not blow up 
below the GUT scale in the NMSSM supplemented by three 5O(10) 10-plets as in Tabled 
The slightly larger value of A here allows extra matter to be at or above the TeV scale. 
This set of benchmarks is displayed in Table 

• NMSSM-with-extra-matter and Lightest Higgs being SM-like near 125 GeV 

This is achieved with A = 0.7 and k = 0.20 — 0.25 such that A does not blow up below 
the GUT scale in the NMSSM supplemented by three 50(10) 10-plets as in Tabled The 
larger values of A and k here requires the extra matter to be close to the electroweak scale. 
This set of benchmarks is displayed in Table [9l 

Note that the input values for A and n are taken at the scale 1 TeV. Their corresponding 
values at Mz are lower so that we are well within the limits given in Tables [JJEJ Within each 
set we have selected three benchmark points which are chosen to illustrate key features of the 
NMSSM Higgs boson which could be used to resolve it from the SM Higgs state in future LHC 
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Table 5: The NLO and NNLO cross-sections for gluon fusion into a SM Higgs boson and the total cross-section 
as defined in Eq. (|4.30[) including an NLO and NNLO gluon fusion cross-section, respectively, for various values 
of the SM Higgs boson mass. 



M H [GeV] 


< L ° [Pb] 


NNLO r_v.l 


[Pb] 


NNLO r„Kl 


123.5 


13.08 


15.44 


0.134 


0.155 


123.6 


13.06 


15.41 


0.135 


0.156 


123.7 


13.04 


15.39 


0.136 


0.157 


123.8 


13.02 


15.36 


0.137 


0.158 


124 


12.97 


15.31 


0.138 


0.160 


124.5 


12.87 


15.18 


0.143 


0.165 


124.6 


12.85 


15.15 


0.144 


0.166 


125 


12.76 


15.05 


0.147 


0.170 


125.8 


12.59 


14.85 


0.154 


0.178 


126 


12.55 


14.80 


0.156 


0.180 


126.2 


12.51 


14.75 


0.158 


0.182 


126.5 


12.45 


14.68 


0.161 


0.185 



searches. This makes 12 benchmark points in total which we label as NMP1 to NMP12. In the 
following we discuss the key features of each of the points in detail under the above four headings. 
Since we focus on large values of A > 0.55 and moderate values of tan (3 = 2, 3 the requirement 
of the validity of perturbation theory up to the GUT scale constrains the low energy value of the 
coupling to be k < 0.3. As a result the heaviest CP-even, CP-odd and charged Higgs states are 
almost degenerate and substantially heavier than the two lightest Higgs scalars and the lightest 
Higgs pseudoscalar for all the benchmark points, as discussed in p3]. Also note that all the 
benchmark points satisfy the WMAP relic abundance. Depending on the decomposition of the 
lightest neutralino and its mass value the main annihilation channels are XiXi ~~ * W + W~ or 
XiXi ~~ QQ (q = b, s,d,c,u). For very light stops t\ with mass near the x? mass, as in the 
benchmark point NMP8, the dominant channel is the coannihilation channel, ~~ > W + b. 

5.1 NMSSM with Lightest Higgs being SM-like 

Table [6] shows three points NMP1, NMP2, NMP3 which satisfy NMSSM perturbativity up to 
the GUT scale and lead to the SM-like Higgs boson always being the lightest one H\ near 
125 GeV. The mass of the second lightest Higgs boson H2 ranges between 129 and 155 GeV. 
Being mostly singlet-like its couplings to SM particles are suppressed enough not to represent 
any danger with respect to the LHC Higgs exclusion limits in this mass range. The benchmark 
points are ordered in terms of increasing stop mixing Xt scaled by the geometric mean of 
the soft SUSY breaking stop mass parameters. The gluino mass is about 700 GeV for all the 
points and recall that the squark and slepton masses not shown are all about 1 TeV. 

NMP1 with (A, k) = (0.64,0.25) and tan/3 = 3 is an example of a benchmark point where a 
SM-like Higgs boson with mass 124.5 GeV can be achieved in the pure NMSSM with relatively 
small mixing, Xt/mf = 1.74. Due to a slightly larger coupling of the SM-like Higgs boson to 
down-type quarks in this scenario, the decay into bb is enhanced compared to the SM and also 
the corresponding branching ratio (normalised to the SM) of R b i = 1.08. As the decay into bb 
contributes dominantly to the total width, its enhanced value leads to smaller branching ratios 
into ZZ,WW, Ryy = 0.94 (V = W,Z). The partial width into photons, however, is enhanced 
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compared to the SM due to the additional SUSY particle loops contributing to the Higgs-77 
coupling with the main enhancement induced by the chargino loops, so that in the end the 
branching ratio into photons is slightly larger than in the SM, i? 77 = 1.03. As for the gluon 
fusion cross-section, which is dominated by the (s)top quark loops at small tan /3, it is almost the 
same as in the SM, with i? CT = 0.97. This is because the Hi couples with SM-strength to the 
top quarks and the effect of the additional squark loops is small. Due to the smaller branching 
ratios into VV, in particular the dominating one into WW, the total Higgs production cross- 
section is estimated to be smaller than in the SM, R atot = 0.92, which is compatible with LHC 
constraints. The heavier stop mass here of 782 GeV is a little uncomfortable from the point of 
view of fine-tuning, but not too bad. 

NMP2 with (A, k) = (0.60,0.18) and tan f3 = 2 shows a 126.5 GeV Higgs which has almost 
SM-like branching ratios into massive vector bosons, Rvv = 1-02. While it will be difficult to 
distinguish this Higgs boson from the SM through the decays into WW, ZZ, the branching ratio 
into 77 is significantly larger than the SM value, -R 77 = 1.20, leading to an enhanced number of 
events in the photon final state, i? (7inci -R 77 = 1.15. The reason for the increased 77 branching 
ratio is an enhanced photonic decay width due to positive contributions from the SUSY loops 
not present in the SM, in particular from the chargino loops, and a reduced total width of 
Rr tot = 0.96. The slight deviation of the branching ratio into bb from the SM value, R^ = 1.05, 
is within our theoretical error estimate. The gluon fusion production cross-section ratio ends up 
being close to the SM value for the same reasons as for NMP1, R agg = 0.96. Due to the SM-like 
decays into massive vector bosons the total Higgs cross-section times branching ratios, which 
is mainly influenced by the branching ratio into WW, is estimated to be very close to the SM 
value, R,j tot = 0.99. The heavier stop mass here of 838 GeV is becoming more uncomfortable 
from the point of view of fine-tuning. 

NMP3 with (A,k) = (0.57,0.20) and tan/3 = 2 shows a 124.6 GeV Higgs boson. Its distin- 
guishing features are a large stop mixing Xt/mj = 2.26 and a considerably enhanced branching 
ratio (normalised to the SM) into photons, i? 77 = 1.42. The large mixing helps to increase the 
Higgs boson mass as in Eq. (|2.8|) and consequently allows the stop masses to be reduced, with 
the heavier stop mass of 686 GeV being apparently more acceptable from the point of view of 
fine-tuning, but at the price of an increased stop mixing angle. The coupling of the SM-like 
Higgs boson to down-type fermions is reduced by about 15% compared to the SM so that the 
total Higgs width dominated by the decay into bb is reduced, Rr tot = 0.78 (normalised to the 
SM). Also the couplings to massive gauge bosons are reduced by 6% implying smaller decay 
rates into WW, ZZ. The branching ratios (normalised to the SM) into WW, ZZ are, however, 
enhanced, Ryy=\.Yl. This is due to the reduced total width. The large enhancement of i? 77 is a 
combination of positive contributions from the SUSY particle loops, mainly charginos, together 
with a reduced total width@ While the coupling to the top quarks is almost SM-like negative 
contributions from the stop quark loops with relatively light stops and a large Higgs coupling to 
these states in this scenario decrease the gluon fusion cross-sectiord so that it amounts only to 
about 77% of the SM value, leading to only ~ 10% more events in the 77 final state compared to 
the SM, Ra incl -R 77 = 1.11, whereas the number of WW events is suppressed, R aincl Rw = 0.88. 
Thanks to the enhanced branching ratios into 77, ZZ and WW, a total cross-section times 
branching ratios still compatible with the LHC searches is achieved, R atot = 0.89. 

4 The phenomenon of an enhanced decay width into photons due to a suppressed decay into bb has also been 
discussed in |101| for the lighter Higgs Hi and in 6 for H2 as well as in [1]. 

5 Note that (s)bottom loops only play a minor role in all our scenarios as we have small values of tan (3. 
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Table 6: NMSSM benchmark points with the lightest Higgs boson being SM-like near 125 GeV. 



Point 


NMP1 


NMP2 


NMP3 


tan 8 


3 


2 


2 


u„ff GeVl 


200 


200 


200 


A 


0.64 


0.6 


0.57 


K 


0.25 


0.18 


0.2 


A x [GeV] 


560 


405 


395 


A K [GeV] 


-10 


-10 


-80 




650 


700 


530 


M- tn [GeV] 


650 


700 


530 


Mi [GeV] 


106 


91 


115 


M 2 [GeV] 


200 


200 


200 


Ms [GeV] 


600 


600 


600 


SM-like Higgs boson 


M Hl [GeV] 


124.5 


126.5 


124.6 




1.03 


1.20 


1.42 


Rvv(Hi) 


0.94 


1.02 


1.12 


R b l{Hi) 


1.08 


1.05 


1.01 


Rr tot (Hi) 


1.05 


0.96 


0.78 




0.97 


0.96 


0.77 


Rincl (-^"l ) 


1.00 


1.15 


1.11 


RinclRw{H\) 


0.91 


0.98 


0.88 


Ra tot (Hl) 


0.92 


0.99 


0.89 


Remaining Higgs spectrum 


M„ 2 [GeV] 


155 


132 


129 


M„ 3 [GeV] 


637 


465 


456 


M Al [GeV] 


128 


116 


168 


Ma 2 [GeV] 


634 


463 


454 


M H ± [GeV] 


626 


454 


447 


Sparticle masses and stop mixing 


m- g [GeV] 


700 


701 


696 


m.± [GeV] 


137 


131 


131 


m.± [GeV] 


281 


284 


284 


m i? [GeV] 


78 


68 


85 


m^o [GeV] 


134 


127 


140 


m^o [GeV] 


201 


169 


178 


m^o [GeV] 


-231 


-232 


-227 


m i? [GeV] 


292 


290 


290 


m 5l [GeV] 


667 


715 


538 


m g , [GeV] 


1014 


1015 


1011 


m fl [GeV] 


548 


587 


358 


m 4 - 2 [GeV] 


782 


838 


686 


Xt/mi 


1.74 


1.86 


2.26 


Relic density 




0.9819 


0.1170 


0.1100 
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5.2 NMSSM with Second Lightest Higgs being SM-like 



Table H shows three points NMP4, NMP5, NMP6 which satisfy NMSSM perturbativity up to 
the GUT scale and lead to the SM-like Higgs boson always being the second lightest one H2 
near 125 GeV. The mass of the lightest Higgs boson H\ is of about 90-96 GeV. We have verified 
in each case that H\ is not in conflict with LEP limits due to its singlet character and hence 
strongly reduced couplings to the SM particles. The gluino mass for these points ranges from 
700 to 875 GeVrJ and, as in all cases, the squark and slepton masses not shown are all about 
1 TeV. 

NMP4 with (A, k) = (0.67,0.10) and tan/3 = 3 shows a 123.8 GeV Higgs boson with almost 
SM-like branching ratios into massive vector bosons, Ryy = 1.01. It differs from the SM Higgs 
by its slightly larger value of branching ratio (normalised to the SM) into photons, -R 77 = 1.06, 
which is due to the positive SUSY particle loop contributions, mainly from chargino loops. The 
combined effect of the sbottom and stop loops in addition to the heavy quark loops is such that 
the gluon fusion production cross-section is SM-like. The total production cross-section times 
branching ratios is near the SM value, R atot = 1.01. With Xt/mi = 1.56, the heavier stop mass 
here of 724 GeV is a little uncomfortable from the point of view of fine-tuning, but not too bad. 
While this Higgs boson will be difficult to be distinguished from the SM Higgs state, the main 
new interesting feature is the appearance of two light bosons with Mu x = 90 GeV and Ma 1 = 86 
GeV which are predominantly composed of singlet S and therefore would have escaped detection 
at LEP. 

NMP5 with (A, k) = (0.66,0.12) and tan /3 = 3 shows a 126.5 GeV Higgs boson which is also 
hard to distinguish from the SM Higgs due to its gauge boson branching ratios being close to 
the SM values, Ryy = 1.06. However the one into 77 is significantly larger, i? 77 = 1.15, which 
is the combined effect of an enhancement due to SUSY particle loops and a smaller total decay 
width, Rr tot = 0.93, mainly due to a smaller decay width into bb as consequence of a reduced 
Higgs coupling to down-type fermions. With the additional suppression of other decay channels, 
in the end the branching ratio into bb is slightly larger than in the SM, R b r = 1.03. The gluon 
fusion production cross-section is only slightly smaller than in the SM, R(j gg = 0.96. Together 
with the enhanced branching ratios into gauge bosons, we estimate the total production Higgs 
cross-section times branching ratios to be close to the SM value, R atot = 1.03, while the number 
of events in the 77 final state is enhanced compared to the SM, R<r incl R~f^ = 1.11. The heavier 
stop mass here of 741 GeV is a little uncomfortable from the point of view of fine-tuning, but 
not too bad. There are two light bosons with Mh 1 = 96 GeV and Ma 1 = 93 GeV which are 
predominantly composed of singlet S and therefore would have escaped detection at LEP. 

NMP6 with (A, k) = (0.55,0.31) and a lower tan/3 = 2 shows a 124.5 GeV Higgs boson which 
can be distinguished from the SM Higgs state by its branching ratios (normalised to the SM) 
into massive gauge bosons of Ryy = 1.10 and in particular by into photons. With i? 77 = 1.39 
the 77 branching ratio shows a very significant enhancement. Once again this is the combined 
effect of a reduction of the total width, -Rrt ot = 0.81, and a larger Higgs to photon coupling due 
to SUSY particle loops. The former results from a reduced Higgs coupling to down-type quarks, 
suppressing the partial width into bb. Nevertheless, the branching ratio for this final state is 
SM-like, R bb = 1.00, due to the additional suppression of other decay channels contributing 

6 The higher value of the gluino mass for the benchmark point NMP6 results from the gaugino mass chosen to 
be M3 = 800 GeV in order to fulfill approximate GUT relations between the soft SUSY breaking gaugino masses. 
This condition is not necessary, however, and choosing a lower value of M3 and hence m§ does not change our 
results in the Higgs sector. 
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Table 7: NMSSM benchmark points with the second lightest Higgs boson being SM-like near 125 GeV. 



Point 


NMP4 


NMP5 


NMP6 


tan 8 


3 


3 


2 


u„ff GeVl 


200 


200 


140 


A 


0.67 


0.66 


0.55 


K 


0.1 


0.12 


0.31 


A x [GeV] 


650 


650 


210 


A K [GeV] 


-10 


-10 


-210 




600 


600 


800 


M- tn [GeV] 


600 


600 


600 


Mi [GeV] 


200 


200 


145 


M 2 [GeV] 


400 


400 


300 


Ms [GeV] 


600 


600 


800 


SM-like Higgs boson 


M„ 2 [GeV] 


123.8 


126.5 


124.5 


Rjj(H 2 ) 


1.06 


1.15 


1.39 


Rvv(H 2 ) 


1.01 


1.06 


1.10 


R b l{H 2 ) 


1.05 


1.03 


1.00 


Rr tot (H 2 ) 


0.99 


0.93 


0.81 




1.00 


0.96 


0.91 


RinclRy/ {H 2 ) 


1.06 


1.11 


1.26 


RinclRw(H 2 ) 


1.01 


1.03 


1.00 


R* tot (H 2 ) 


1.01 


1.03 


1.02 


Remaining Higgs spectrum 


M Hl [GeV] 


90 


96 


90 


M„ 3 [GeV] 


654 


656 


325 


M Al [GeV] 


86 


93 


249 


Ma 2 [GeV] 


655 


656 


317 


M H ± [GeV] 


643 


645 


312 


Sparticle masses and stop mixing 


m- g [GeV] 


699 


699 


875 


m.± [GeV] 


182 


182 


114 


m.± [GeV] 


438 


438 


342 


m i? [GeV] 


69 


78 


80 


m^o [GeV] 


173 


175 


163 


m^o [GeV] 


233 


238 


-169 


m^o [GeV] 


-241 


-239 


197 


m i? [GeV] 


438 


438 


343 


m 5l [GeV] 


619 


617 


822 


m u [GeV] 


1013 


1013 


827 


m fl [GeV] 


517 


483 


549 


m h [GeV] 


724 


741 


892 


Xt/mi 


1.56 


1.89 


-1.83 


Relic density 


£lh 2 


0.0999 


0.1352 


0.1258 
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to T t ot- In this scenario, the squark loops reduce the gluon fusion cross-section by about 9%, 
R agg = 0.91, resulting in an estimate of the total cross-section times branching ratios which due 
to the enhanced branching ratios into vector bosons is very close to the SM, R atot = 1-02. The 
number of photon final state events, however, is enhanced by considerable 26% compared to the 
SM. The mixing of Xt/m^ = —1.83 leads to stop masses of 549 and 892 GeV, the latter making 
the scenario more uncomfortable from the point of view of fine-tuning. Interestingly, in this 
case there is only one light boson with Mh 1 = 90 GeV, which is predominantly singlet-like and 
therefore would have escaped detection at LEP. 

5.3 NMSSM-with-extra-matter and Second Higgs being SM-like 

Table El shows three points NMP7, NMP8, NMP9 which satisfy perturbativity up to the GUT 
scale providing that the NMSSM is supplemented by three S'O(IO) 10-plets as in Table [3l The 
only slightly larger value of A = 0.68 — 0.69 allows the extra matter to be at or above the TeV 
scale. It therefore may play no role in any LHC phenomenology, apart from the indirect effect 
of allowing A to be a bit larger than allowed by the usual requirement of perturbativity in the 
NMSSM with no extra matter. All these points lead to the SM-like Higgs boson always being 
the second lightest one H2 near 125 GeV. The points are ordered in terms of increasing stop 
mixing AQ scaled by mf. In two of the benchmark scenarios H\ is light with masses of 83 and 
93 GeV, respectively, in one scenario it is almost degenerate with H2 with a mass of 124 GeV. 
In all these cases we have verified that such light bosons are sufficiently singlet-like that they 
would have been produced with sufficiently low rates in order not to be in conflict with LEP 
and LHC searches. The gluino mass ranges between about 700 and 880 GeV. Recall that the 
squark and slepton masses not shown are all about 1 TeV. 

NMP7 with (A, «) = (0.68, 0.06) and tan/3 = 2 involves a 124.5 GeV Higgs boson which will be 
difficult to be distinguished from the SM Higgs by branching ratios into massive vector bosons 
(normalised to the SM) of Rvv = 1-06. The reduced coupling to down-type quarks leads to 
a smaller decay width into bb and hence a smaller total width than in the SM, Rr tot = 0.90. 
The combined effect of a reduced decay into bb and a positive contribution of, in particular, 
chargino loops to the Higgs coupling to photons leads to a significantly enhanced branching 
ratio into photons compared to the SM with i? 77 = 1.17, which can be used to distinguish this 
NMSSM scenario from the SM case through the enhanced number of events in the photon final 
state, R amcl Rjj = 1.20. The H2 coupling to up- type quarks is SM-like so that the gluon fusion 
production dominated by (s)top quark loops is not altered. With a small positive contribution 
from stop loops this leads to a 3% enhancement, Ra aa = 1-03, and implies an estimated total 
Higgs production cross-section times branching ratio (normalised to the SM value) of R atot = 
1.10. The stop masses are almost 800 GeV making this point appear to be somewhat fine-tuned. 
However, the absence of any stop mixing AQ/mj = 0, which is the main feature of this point, 
serves to ameliorate the fine-tuning. Note that, in the MSSM, a 125 GeV Higgs with zero stop 
mixing would be impossible. In addition there are two light bosons with Mh x = 93 GeV and 
M J 4 1 = 99 GeV which are predominantly singlet-like and therefore would have escaped detection 
at LEP. 

NMP8 with (A, n) = (0.69,0.125) and tan/3 = 2 involves a 126.2 GeV Higgs boson. The 
main features of this scenario are reduced branching ratios not only into massive gauge bosons, 
Rvv = 0.87, but also into photons, i? 77 = 0.87, and particularly light stop masses with the 
lighter = 104 GeV. The branching ratio reductions are due to an enhanced decay width into 
bb because of a slightly larger coupling to down-type quarks than in the SM. In addition, the 
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Table 8: NMSSM-with-extra-matter benchmark points with second lightest Higgs being SM-like near 125 GeV. 



Point 


NMP7 


NMP8 


NMP9 


tan 8 


2 


2 


3 


u„ff fGeVl 


200 


190 


200 


A 


0.68 


0.69 


0.68 


K 


0.06 


0.125 


0.2 


A x [GeV] 


500 


420 


600 


A K [GeV] 


-100 


-100 


-10 


M & [GeV] 


750 


400 


500 


M- tR [GeV] 


750 


230 


500 


Mi (GeV) 


200 


120 


135 


M 2 (GeV) 


400 


335 


200 


M :i (GeV) 


800 


800 


600 


SM-like Higgs boson 


M„ 2 [GeV] 


124.5 


126.2 


125.8 


Rjj(H 2 ) 


1.17 


0.87 


1.78 


Rvv(H 2 ) 


1.06 


0.87 


1.41 


R b l{H 2 ) 


1.02 


1.12 


0.84 


Rr tot (H 2 ) 


0.90 


1.06 


0.53 




1.03 


1.01 


0.79 


RinclRy/ {H 2 ) 


1.20 


0.86 


1.39 


RinclRw(H 2 ) 


1.09 


0.87 


1.11 


R* tot (H 2 ) 


1.10 


0.87 


1.13 


Remaining Higgs spectrum 


M Hl [GeV] 


93 


83 


124 


M„ 3 [GeV] 


495 


434 


641 


M Al [GeV] 


99 


139 


119 


Ma 2 [GeV] 


503 


438 


640 


M H ± [GeV] 


485 


422 


629 


Sparticle masses and stop mixing 


m- g [GeV] 


882 


877 


695 


m.± [GeV] 


179 


159 


137 


m.± [GeV] 


439 


379 


281 


m i? [GeV] 


66 


76 


81 


m^o [GeV] 


160 


120 


143 


m^o [GeV] 


231 


191 


185 


m^o [GeV] 


-248 


-233 


-237 


m i? [GeV] 


440 


380 


291 


m 5l [GeV] 


779 


379 


514 


m g , [GeV] 


1022 


999 


1011 


m h [GeV] 


786 


104 


391 


m 4 - 2 [GeV] 


797 


442 


634 


Xt/mi 





1.48 


1.77 


Relic density 


n/i 2 


0.0973 


0.1101 


0.1031 
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decay into photons is reduced because of smaller H2 couplings to up-type quarks and negative 
stop quark loop contributions in this scenario. In the gluon fusion production the positive 
contributions from the stop quark loops counterbalance the reduction of the cross-section because 
of the smaller H2 coupling to top quarks, so that the ratio to the SM cross-section is nearly 
one, R(j gg = 1.01. However the smaller branching ratios into WW,ZZ and 77 drive the total 
cross-section times branching ratios down to R atot = 0.87 (normalised to the SM). This is still 
compatible with the present LHC searches though, according to our criteria. With X-t/m$ = 1.48, 
the interesting feature of this scenario is the very light stop which is close enough in mass to 
the lightest neutralino in order to dominantly decay via the flavour changing neutral current 
decay into charm and Xi- Such a light stop is not excluded in this decay channel by the present 
LHC searches. The correct amount of dark matter relic density is here achieved through t\x\ 
co-annihilation into W + b. Also in this scenario we have two light bosons with Mh 1 = 83 GeV 
and M J 4 1 = 139 GeV which are predominantly singlet-like and hence not in conflict with LEP 
and LHC exclusion limits. 

NMP9 with (A,k) = (0.68,0.20) and tan/3 = 3 involves a 125.8 GeV Higgs boson. The main 
features of this scenario are the near degeneracy of the physical SM-like Higgs with a lighter Higgs 
boson of mass Mh 1 = 124 GeV as well as considerably enhanced branching ratios (normalised 
to the SM) into photons, _R 77 = 1.78, and into massive gauge bosons, Rvv = 1.41. Due 
to the near degeneracy of H% and H2, and the associated large mixing, each state contains a 
significant singlet component, with the most SM-like Higgs H2 coupling only with about 60% 
of the SM strength to down-type quarks so that both the decay width into bb and the total 
width are reduced, R b i = 0.84 and Rr tot = 0.53, and hence the branching ratios into gauge 
bosons enhanced. As the H2 coupling to top quarks is also reduced the gluon fusion production 
cross-section goes down supplemented by a small negative contribution from stop loops, so that 
in the end we have R<j gg = 0.79. Nevertheless, due to the enhanced branching ratios, the 
combined cross-section times branching ratios is compatible with LHC searches as R atot = 1.13. 
Particularly interesting is the still considerably enhanced number of photon final state events 
compared to the SM, R aincl R~~ = 1.39, which is a distinctive feature of this scenario. Also 
the Hi state is consistent with LHC searches as its couplings to up-type quarks are reduced so 
that it cannot be produced with sufficiently large rate to yield a second signal at the LHC. In 
addition there is a singlet-dominated boson of mass Ma x = 119 GeV. With Xt/M m - = 1.77, 
and the heavier stop mass here of 634 GeV, and a gluino mass of 695 GeV, there is relatively 
little fine-tuning for this point. 

5.4 NMSSM-with-extra-matter and Lightest Higgs being SM-like 

Table M shows three points NMP10, NMP11, NMP12 which satisfy perturbativity up to the 
GUT scale providing that the NMSSM is supplemented by three 5O(10) 10-plets as in Table [2 
The larger value of A = 0.7, held fixed for all these points, requires the extra matter to be close 
to the electroweak scale. It may therefore be expected to play a role in LHC phenomenology. 
All these points lead to the SM-like Higgs always being the lightest one, Hi, near 125 GeV. 
The points are ordered in terms of increasing stop mixing Xt scaled by mj. The second lightest 
Higgs boson has masses in the range of about 130 to 145 GeV and, due to its singlet nature, 
we have checked that it is in accordance with present LHC exclusion limits. The gluino mass 
for these benchmark points is about 700 GeV and recall that the squark and slepton masses not 
shown are all about 1 TeV. 

NMP10 with (A, k) = (0.70,0.20) and tan/3 = 2 involves a 123.6 GeV Higgs which is SM-like 
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in its decays into massive vector bosons and b quarks with the normalised branching ratios given 
by Rvv = 1-01 and R b i = 1.04. The branching ratio into photons, however, is enhanced with 
(normalised to the SM) i? 77 = 1.15. This is a consequence of positive contributions from chargino 
loops to the fli77 coupling. The cross-section ratio R<j gg = 1-07 is somewhat larger than in the 
SM because of positive contributions from the stop loops. Together with the SM-like branching 
ratios into massive gauge bosons and an enhanced one into 77, in the end the estimated total 
Higgs cross-section times branching ratios is not far from the SM value, R Utot = 1-08. The 
number of photon events, however, is enhanced by significant 21% compared to the SM. The 
stop masses are 630 and 645 GeV with the gluino mass 686 GeV, which does not require fine- 
tuning, since there is no stop mixing Xt/m^ = 0. This of course would be impossible in the 
MSSM. The reason it is possible here is the large A = 0.7, and low tan/3 = 2, which together 
provide a nice tree- level Higgs mass contribution as in Eq. (|1.2p . The heavier Higgs H2 is not 
too far away at 133 GeV but, of course, is singlet dominated and hence hard to discover. 

NMP11 with (A,k) = (0.70,0.20) and tan/3 = 2 involves a 123.5 GeV Higgs. Here we have 
a scenario with a strong singlet-doublet mixing, this time for the lightest Higgs boson H\. As 
a consequence, its couplings to down-type quarks are reduced implying a smaller decay width 
into bb leading to a total width of only about 40% of the SM value, Rr tot = 0.40. The branching 
ratios into massive gauge bosons are therefore significantly enhanced, Rvv = 1-48, despite 
Hi also having reduced couplings to gauge bosons. The most interesting feature, however, is 
the enhancement of the branching ratio into photons by a factor 2, i? 77 = 2.00, due to the 
combined effect of the smaller total width and the positive chargino loop contributions which 
counterbalance the loss in the Hi coupling to 77 because of a reduced coupling to top quarks. 
The cross-section ratio R<j gg = 0.71 results from the reduced top loop contributions which cannot 
be made up for by positive contributions from the squark loops. While the number of WW final 
states and the total cross-section times branching ratios, dominated by the decay into WW, end 
up near their corresponding SM values, the number of photonic events is enhanced by as much 
as 39%, providing a distinctive feature for this benchmark point. The nearby H2 state with 136 
GeV represents no danger with respect to LHC exclusion limits as its couplings to SM particles 
are reduced implying small enough production rates to be safe. With X t /m^ = 1.33, and a 
heavier stop mass of 620 GeV, and a gluino mass of 696 GeV, this is a natural point without 
much fine-tuning. 

NMP12 with (A, k) = (0.70,0.25) and tan /3 = 3 involves a 123.7 GeV Higgs boson which 
is harder to distinguish from the SM Higgs by branching ratios into massive vector bosons 
(normalised to the SM) of Rvv = 1-09. However with i? 77 = 1.24 the branching ratio into 
photons is significantly different from the SM and enhanced once again because of positive 
chargino loop contributions. As the Hi coupling to top quarks is the same as in the SM and 
the marginal contributions from squark loops add up to almost zero, the gluon production 
cross-section is the same as in the SM, R agg = 1-00. With the enhanced branching ratios into 
vector bosons the total production times branching ratios ends up to be larger than in the SM, 
Ratot = 1-10, and the number of photon final state events is enhanced by significant 23%. The 
normalized mixing Xt/m^ = 1.57 features stop masses of 413 and 624 GeV. With the gluino 
mass of 696 GeV these are quite natural values which does not require fine-tuning. 
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Table 9: NMSSM-with-extra-matter benchmark points with the lightest Higgs being SM-like near 125 GeV. 



Point 


NMP10 


NMP11 


NMP12 


tan /3 


2 


2 


3 


Ur-n- fGeVl 


200 


200 


200 


A 


0.7 


0.7 


0.7 


K 


0.2 


0.2 


0.25 


A x [GeV] 


405 


405 


560 


A K [GeV] 


-10 


-10 


-10 




600 


510 


500 


M in [GeV] 


600 


510 


500 


M 1 [GeV] 


120 


120 


125 


M 2 [GeV] 


200 


200 


200 


M 3 [GeV] 


600 


600 


600 


SM-like Higgs boson 


M Hl [GeV] 


123.6 


123.5 


123.7 


R~ n (H\) 


1.15 


2.00 


1.24 


Rvv{Hi) 


1.01 


1.48 


1.09 


Rbb(Hi) 


1.04 


0.82 


1.01 


Rr tot {Hi) 


0.98 


0.40 


0.90 




1.07 


0.71 


1.00 


Rind R^fy (-^1 ) 


1.21 


1.39 


1.23 


Rincl Rvv(Hl) 


1.07 


1.03 


1.09 


R* tot (Hi) 


1.08 


1.06 


1.10 


Remaining Higgs spectrum 


M„ 2 [GeV] 


133 


136 


144 


M H3 [GeV] 


468 


458 


629 


M Al [GeV] 


129 


129 


133 


Ma 2 [GeV] 


469 


459 


626 


M H ± [GeV] 


456 


445 


615 


Sparticle masses and stop mixing 


nig [GeV] 


686 


696 


696 


m.± [GeV] 


131 


131 


137 


m-± [GeV] 


284 


283 


281 


m^o [GeV] 


82 


82 


83 


m-o [GeV] 


137 


137 


141 


m^o [GeV] 


175 


175 


198 


m^o [GeV] 


-241 


-240 


-236 


m^o [GeV] 


291 


291 


293 


m 5l [GeV] 


622 


525 


516 


m g , [GeV] 


1013 


1011 


1011 


m fl [GeV] 


630 


445 


413 


m t - 2 [GeV] 


645 


620 


624 


X t /m f 





1.33 


1.57 


Relic density 


fi/i 2 


0.1135 


0.1161 


0.1050 
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6 Summary and Conclusion 



The recent ATLAS and CMS indications of a SM-like Higgs boson near 125 GeV are consistent 
not only with the SM Higgs but also a SUSY Higgs. If it is a SUSY Higgs, the question is what 
sort of SUSY theory is responsible for it, and what part of parameter space of the SUSY theory 
does it correspond to? In this paper we have been guided by both naturalness and minimality 
to consider the possibility that the SM-like Higgs arises from the NMSSM. Minimality alone 
would suggest the MSSM, however a 125 GeV Higgs can only be achieved in the MSSM with 
stops necessarily having some mixing and mass parameters above 1 TeV, which leads to some 
considerable amount of tuning. By contrast, the NMSSM, which includes only one additional 
singlet superfield S, allows a 125 GeV Higgs to arise with all stop masses and mixing below 1 TeV. 
Thus the NMSSM appears to be the best compromise between naturalness and minimality that 
can account for a 125 GeV SM-like Higgs boson, with the SM itself excluded on naturalness 
grounds. 

Given that the NMSSM is a well-motivated SUSY theory, the next question is what part of 
NMSSM parameter space does the SM-like Higgs boson near 125 GeV correspond to? Clearly 
on naturalness grounds we are led to consider large values of A and low values of tan /?, in order 
to maximise the tree-level contribution to the Higgs mass in Eq. (|1.2|) . However A cannot be 
too large otherwise it would blow up below the GUT scale which we regard as unacceptable. 
Therefore it is important to know just how large A can be, as a function of k and tan f3. We 
have studied this question using two-loop RGEs and found that it cannot be as large as is often 
assumed in the literature. However, if the SUSY desert contains extra matter, then A can be 
increased, depending on the precise mass threshold of the extra matter, which is a new result of 
this paper. 

Assuming the NMSSM with large values of A and low values of tan j3, which may or may 
not require extra matter, another key requirement for us is having stop masses and mixings 
below 1 TeV, on naturalness grounds. Since the LHC has already placed strong limits on the 
first and second family squarks we are therefore led to consider a general NMSSM, i.e. not the 
constrained NMSSM. Given the rich parameter space of the NMSSM, it is practically impossible 
to perform scans over all NMSSM parameter space to identify which regions are consistent with a 
125 GeV Higgs. Instead we have adopted a different strategy in this paper, namely that of using 
benchmark points. Each benchmark point in the NMSSM parameter space allows everything 
to be calculated to the standard of precision that is currently available. For each benchmark 
point, the nature and properties of the SM-like NMSSM Higgs boson can be scrutinised and 
compared to the SM expectation, which could allow the NMSSM Higgs to be resolved from the 
SM Higgs boson by future LHC searches and measurements. It is possible for experimentalists 
to use these benchmark points to test the NMSSM in different regions of the parameter space 
by performing studies beyond what, as mere theorists, we are capable of. 

Motivated by these considerations we have proposed four sets of benchmark points corre- 
sponding to the SM-like Higgs being the lightest or the second lightest in the NMSSM or the 
NMSSM-with-extra-matter. Points NMP1, NMP2, NMP3 are for the NMSSM with the lightest 
Higgs being SM-like near 125 GeV. Points NMP4, NMP5, NMP6 are for the NMSSM with the 
second lightest Higgs being SM-like near 125 GeV. Points NMP7, NMP8, NMP9 are for the 
NMSSM-with-extra-matter at or above the TeV scale where the second lightest Higgs is SM-like 
near 125 GeV. Points NMP10, NMP11, NMP12 are for the NMSSM-with-extra-matter close to 
the electroweak scale and the lightest Higgs being SM-like near 125 GeV. For each of these points 
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we have discussed how the NMSSM Higgs boson near 125 GeV may be distinguished from the 
SM Higgs state in future LHC searches. We have considered both the gluon gluon cross-section 
and the branching ratios into WW, ZZ and two photons. We also considered a weighted total 
Higgs cross-section times branching ratios as a measure of how many Higgs events would be 
observed in current searches. 

Each benchmark point, with different values of Higgs cross-sections and branching ratios, 
acts like a little beacon which in principle can allow experiment to guide us to the correct part 
of the NMSSM parameter space. It is worth emphasising, however, that the values of Rvv and 
R b l quoted here are only accurate to a few percent, as they are obtained by using two different 
programs for the calculation of the branching ratios, NMSSMTools for the NMSSM and HDECAY 
for the SM case, which use different approximations in the calculation of the various partial 
widths. The differences are always less than or equal to 5% for the quoted values, while the 
maximum deviations for i? 77 differ by up to 10% . Although these theoretical uncertainties 
are well within the expected experimental error for the measurement of these branching ratios 
at the LHC, it does mean that resolving the NMSSM Higgs from the SM Higgs is going to be 
challenging. However, in some cases the Higgs branching ratio into two photons may easily 
be enhanced by a factor of two in the NMSSM, as seen in NMP9 and NMP11. These very 
large enhancements result in part from the reduction in the bottom quark-antiquark decay rate 
due to the reduced amount of component in the SM-like Higgs boson. Moreover significant 
enhancements in the two photon channel are also seen for other points where the bb decay rate 
is not reduced, and this is typically due to the effect of chargino loops. 

In general we have taken the stop masses and mixing to be as low as possible, the extreme 
example of this being NMP8 where the lightest stop is only 104 GeV. We find it remarkable that 
such light stops are still allowed by current searches, although we would expect stops below 1 
TeV, as assumed here, to be discovered soon at the LHC. We look forward to the forthcoming 
LHC run in 2012 with much excitement. 
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Appendix 



A Renormalisation Group Equations to Second Order 

The running of the gauge couplings from the GUT scale to the EW scale is determined by a set 
of RGEs. In our analysis, we use two-loop RGEs for the gauge and Yukawa couplings. Retaining 
only A, k, ht, h b and h T the two-loop RGEs in the NMSSM can be written as 



dgi 




dt 


(4tt) 2 


d\ 


A 


~dt 


(4vr) 2 


dn 


K 




(4tt) 2 



3 3 {2) i 
4A 2 + 2k 2 + 3h 2 + 3h 2 b +h 2 T - 3g\ - -gj + 

5 (47TJ Z 



6{k 2 + A 2 ) + 



47T 



18A 2 5 2 + ^X 2 g 2 
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dh t 


h t 


dt ~ 


(47T) 


dhf, 


h 


~~dt 


(47T) 


dh T 


h T 


~dt 


(47T) 



A 2 (12A 2 + 24k 2 + 18/i? + 18/i 2 + 6/i 2 ) - 24k 



\2 , cl.2 , 1,2 16 2 o 2 13 2 

A + 6/i t + h b - —g 3 - 3 5 2 " J^9i + 7^2 



\2 1 7,2 1 c t2 . 7,2 16 2 o 2 ^ 2 , ^6 

A + /i t + 6h b + /i T - y g 3 - 3g 2 - — 5l + 



o /?( 2 ) 

A 2 + 3/l 2 +4 ^2_3 5 2__ 5 2 + _^_ 



(A.l) 



where i = \n[Q/Mx\, the index i runs from 1 to 3 and is associated with U(1)y, SU(2)w 

(2) (2) (2) (2) 

and SU(3)c gauge interactions; (3 X , f3 ht , and are the two-loop contributions to the 
corresponding /3-functions. 

(2) (2) (2) (2) 

In the NMSSM the two-loop /3-functions of the gauge couplings, j3 x , j3 ht , (3^ and (3 h J are 
given 

ft 
£2 
ft 



-3 + 



1 



1 + 



16vr 2 
1 



14 5 2 + 9g 2 + H 5 2 - Ah 2 



Ah 2 



16vr 2 



24<? 2 + 25<?| + ^ 



6/i 2 



2h 2 - 2A 2 



(A.2) 



33 1 
~5~ + 16vr 2 



'88 2 27 2 199 2 26 2 14,. 

T 53 + Y 52 + ^ ~ Y ht ~ Y h > 



7 Note that the two-loop RGE results differ somewhat from the original two-loop RGE results first obtained 
in the NMSSM by King and White |13| . which is why we have written them out explicitly here. Of course the 
two-loop RGE results including extra matter are new results presented in this paper for the first time. 
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o(2) 
Ph t 
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207 
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15 4 27 4 
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In the NMSSM with three extra S'O(IO) 10-plets the two-loop /3-functions of the gauge 
couplings, (3^\ ft]®, /3f^ and f3^ can be written as 
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